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The isospin-dependence of the nuclear equation of state (EOS) is one of the most im- 
portant but very poorly known property of neutron-rich matter |T|]. Its determination in 
laboratory-controlled experiments has profound implications to the study of the structure 
and evolution of many astrophysical objects 0]. Nuclear reactions induced by stable neutron- 
rich nuclei and/or radioactive beams provide a means to extract useful information about the 
isospin-dependence of the nuclear EOS and to explore novel phenomena in nuclear matter at 
extreme isospin asymmetries. A number of dedicated experiments to study the isospin de- 
pendence of the nuclear EOS have been performed/planned at several available radioactive 
beam facilities and the future Rare Isotope Accelerator (RIA) . For these experiments to 
be fruitful it is important to first understand well the role of the isospin degree of freedom 
in nuclear reaction dynamics. Moreover, theoretical predictions on the sensitivity of exper- 
imental observables to the isospin dependence of the nuclear EOS are useful for explaining 
available data and planning new experiments. Recently, several useful observables have been 
identified in neutron-rich heavy-ion collisions at intermediate energies. These include the 
isospin fractionation P^TU|, neutron to proton ratio of heavy residues [|n|] or projectile-like 



fragments [0,|T3] and the neutron-proton differential flow Most of these observables 

make use of the relative multiplicities and kinetic energies of mirror nuclei. Information 
about the isospin dependence of the nuclear EOS can also be obtained from studying free 
nucleons as neutrons and protons have the opposite symmetry potentials in nuclear medium. 
In a recent work, one of the present authors has shown that the neutron-proton differential 
flow is a rather useful probe of the isospin dependence of the nuclear EOS. The method 
utilizes constructively both the isospin fractionation and the nuclear collective flow as well 
as their sensitivities to the nuclear EOS. It, however, requires measuring the neutron and 
proton collective flow simultaneously. Experimentally, it is easier to identify charged par- 
ticles and measure their momenta accurately, thus observables using charged particles only 
are more useful. In this work, within an isospin-dependent transport model for nuclear re- 
actions involving neutron-rich nuclei, we explore the sensitivity of proton collective flow to 
the isospin dependence of the nuclear EOS in heavy-ion collisions at intermediate energies. 



We examine both the first-order transverse flow and the second-order differential eUiptic 
flow as functions of transverse momentum. It is found that the differential elliptic flow of 
midrapidity protons, especially at high transverse momenta, is very sensitive to the isospin 
dependence of the nuclear EOS. This sensitivity is much stronger than that found in the 
flrst-order transverse collective flow around the projectile and target rapidity. 

At present, nuclear many-body theories predict vastly different isospin dependence of the 
nuclear EOS depending on both the calculation techniques and the bare two-body and/or 
three-body interactions employed, see e.g., [|14| -[17||. Various theoretical studies (e.g., [[I8i |l9|) 
have shown that the energy per nucleon e(p, 6) in nuclear matter of density p and isospin 
asymmetry parameter 6 deflned as 

S = {Pn- Pp)/{pn + Pp) (1) 

can be approximated very well by a parabolic function 

e(p,5) = e(p,0) + 5(p)-5l (2) 

In the above e(p, 0) is the EOS of isospin symmetric nuclear matter and S{p) is the symmetry 
energy at density p. The EOS should also be momentum-dependent, however, for this 
exploratory study we shall negelect this dependence. We shall use for isospin-symmetric 
nuclear matter a stiff EOS with Kq = 380 MeV which can reproduce the transverse flow 



data equally well as a momentum-dependent soft EOS with Kq = 210 MeV pO| , pT| . This 
choice will not affect our conclusions since the elliptic flow has been shown to be insensitive 
to the momentum-dependence of the EOS in heavy- ion collisions at beam energies below 
about 1.5 GeV/nucleon if Kq = 380 MeV is used ||22|. The form of the symmetry energy 



as a function of density is rather strongly model dependent [^]. Very divergent predictions 
on the isospin dependence of the nuclear EOS by various many-body theories have lead to 
vastly different forms for the S{p). We adopt here a parameterization used by Heiselberg 



and Hjorth- Jensen in their studies on neutron stars ||24 



Sip) = Soipo)-u\ (3) 



where u = p/po is the reduced density and 5*0 (po) is the symmetry energy at normal nuclear 
matter density po- The value of 5*0 (po) is known to be in the range of about 27-36 MeV from 
analyzing atomic masses [p5|-|27| . By fitting the result of variational many-body calculations 



by Akmal et al |]T5[, Heiselberg and Hjorth- Jensen found the values of S{po) =32 MeV and 
7 = 0.6. However, as shown by many other authors previously |^ and more recently by 
Brown ||l6l using other approaches, the extracted value of 7 varies widely, even its sign 



is undetermined. Therefore,in this work we take the 7 as a free parameter and study its 
influence on the experimental observables in neutron-rich heavy-ion collisions. Further, we 
take a constant value of 30 MeV for 6*0 (po) in the following study. 

It is well known that the symmetry energy has a kinetic and a potential contribution 

5(p) = ^e°«i(2i-l) + V^2 (4) 

where e'^ is the Fermi energy in symmetric nuclear matter at normal density and V2 is the 
potential contribution. Within our transport model, the kinetic part of the symmetry energy 
is simulated by using different Fermi momenta for neutrons and protons according to the 
local Thomas-Fermi approximation. While the potential part is taken into account by using 
the symmetry potential. Corresponding to Eq. § the symmetry potential energy density is 

Wasy = V2p6\ (5) 

and the single-particle symmetry potential V^'^y can be obtained from 



V^y = ^TT^ = (Soil - IW + 4.2^3)52 ± ^Sou^ _ 12.7^2/3)5 (6) 

OPq 

where and "-" are for q = neutron and q = proton, respectively. For small isospin 
asymmetries and densities near po the above symmetry potentials reduce to the well-known 
Lane potential which varies linearly with 6 . Shown in Fig. 1 are the symmetry potentials 
at 6=0.2 with the 7 parameter of 0.5, 1.0 and 2.0 for neutrons (upper branch) and protons 
(lower branch), respectively. Generally, the repulsive/attractive symmetry potential for 
neutrons/protons increases with density. The isospin dependence of the nuclear EOS can 
be characterized by the curvature of the symmetry energy at po 



K,ym ^ 9pl^p^\,=,, = 95o(po) ■ 7(7 - 1). (7) 



The slope and curvature of the symmetry potential depends strongly on the parameter 7. We 
consider in the following two cases with 7 = 0.5 and 2, the corresponding Ksym parameter 
is —68 MeV and 540 MeV, respectively. This choice of the parameter 7 allows us to explore 
a large uncertain range of the isospin dependence of the nuclear EOS. 

Our study is based on the isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) 
transport model (e.g., [|,^]). In this model protons and neutrons are initialized in coordi- 



nate space according to their density distributions predicted by the relativistic mean-field 
(RMF) theory [^,^. In momentum space they are initialized in Fermi spheres with in- 
dividual neutron or proton Fermi momentum dertermined by its local density according 
to the Thomas-Fermi approximation. The isospin-dependent reaction dynamics is included 
through isospin-dependent nucleon-nucleon scatterings and Pauli blockings, the symmetry 
potential V^sy and the Coulomb potential for protons. A Skyrme-type parameterization 
is used for the isoscalar potential. For a review of the model, we refer the reader to ref. 

However, a few comments relevant to the present study are necessary here. We did not 
include in the EOS the curvature term which is important for reproducing surface proper- 
ties of static nuclei and keeping them stable. This is because we can not find a universal 
coefficient of this term that can reproduce the neutron and proton density distribtions of 
all nuclei up to mass 200. Instead, we used the predictions by the RMF theory to initialize 
nucleons. It would be ideal if one can use the same model with the same interaction to first 
reproduce perfectly the initial neutron and proton density and momentum distributions, in 
particular, the neutron skins of heavy nuclei, then study the reactions between them. What 
we have done sofar is thus a compromise. Nevertheless, the best available knowledge about 
the inital neutron and proton distributions is used, and we found that the initial distribu- 
tions for stable nuclei are kept stable upto about 150 fm/c with energetic particle emissions 
less than about 1 percent in our approach. Thus the distributions are sufficiently stable over 
a time period that is long enough for studying the collective flow which is mainly generated 



in the early stage of the reaction. 

First, to understand the role of the symmetry potential in the reaction dynamics we have 
studied the central density as a function of time for the reaction of Sn+Sn. To the first order 
of approximation, there is essentially no difference in the evolution of the central density 
with 7 =0.5 and 2 since the symmetry potential is rather small compared to the isoscalar 
nuclear potential. We have also investigated the relative roles of the symmstry and Coulomb 
potentials in the reaction dynamics by turning on or off one of them. We found that the 
Coulomb potential dominates over the attractive symmetry potential for protons for the 
reaction considered here. However, without including the symmetry potential the Coulomb 
potential along will lead artifically to more proton emissions than neutrons in the early 
stage of the reaction. Moreover, stable nuclei can not be kept stable using only the Coulomb 
potential without the symmetry potential in the transport model. We thus either turn on or 
off simultaneously both the Coulomb and symmetry potential in our following studies. To 
evaluate the degree of isospin asymmetry obtained and its correlation with density we show 
in Fig. 2 a scatter plot of the local density p{cell)/pQ and the isospin asymmetry 6{cell) of 
all non-vacant cells of 1 fm^ volume each. The plot is made by using the parameter 7=0.5 
for the Sn+Sn reaction at the time instant of 30 fm/c when the central density is close to 
its maximum. It is seen that the high density region centers around the isospin asymmetry 
of the reaction system of 6q ~ 0.19 as one would expect. Cells with lower densities extends 
in 6 toward the neutron-rich side far away from 6q due to both the neutron-skins of the 
colliding nuclei and the isospin fractionation during the reaction. The latter is an unequal 
partitioning of the neutron to proton ratio N/Z of asymmetric nuclear matter between low 
and high density regions. It is energetically favorable for the asymmetric nuclear matter to 
be separated into a neutron-rich low density phase and an isospin-symmetric high density 
one p2|-p^. Since the symmetry potential is an increasing function of 6, the symmetry 
potentials at densities less than po but with larger isospin asymmetries are also appreciable 
as in the higher density region around the Sq. Although this symmetry potential itself has 
almost no effect on the compression and expansion of the central region, it does have a 
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significant effect on the emission of particles and their kinetic energy spectra. As seen from 
Fig. 1, at densities below about po? the attractive symmetry potential for protons is stronger 
with the parameter 7 =0.5 than 7 = 2.0. As a result, more protons are expected to become 
unbound with 7 =2. Shown in Fig. 3 are the unbound proton rapidity distributions at the 
freeze-out for the Sn+Sn reaction. The unbound protons are identified as those with local 
densities less than po/8 and the freeze-out time is found to be about 100 fm/ c for the reaction 
considered here. Appreciable effects of the symmetry potential on the rapidity distribution of 
free protons are seen over the whole rapidity range. Significantly more particles are emitted 
with 7 = 2., especially around the midrapidity. This observation is consistent with the total 
and differential isospin fractionation discussed in detail in refs. P JT0| , p!3[ . 



We now turn to the exploration of the isospin dependence of the nuclear EOS by studying 
both the first-order transverse collective flow and the second-order differential elliptic flow of 
unbound protons. Collective flow has been well known for its ability to reveal the pressure 



gradient during nuclear reactions, see e.g., p7|-|39|. In the local rest frame of matter the 



collective flow velocity v, energy density e and pressure p are related by the relativistic 
Euler equation 



d 

{e + p)—v = -Vp. (8) 

The pressure gradient Vp plays the role of a driving force for the acceleration of the collective 
flow velocity v. For isospin asymmetric nuclear matter the pressure p has a kinetic part Pkm, 
an isospin-independent interaction part po and an isospin-dependent one Pasy, i-e., 

P = Pkin +P0 +Pasy (9) 

At temperature T the kinetic contribution pkm is 

oo \ 3 



1 oo \ 3 

Rin = Tp{l + bni^r [(1 + 5)^-^" + (1 - Sy^''] }, (10) 

^ n=l ^ 

where the b'^s are the inversion coefficients given in ref. and 



is the thermal wavelength of nucleons with mass m. The contribution from the isospin- 
independent nuclear interaction po is 

Po = ^apo(fr + y^Po(fr\ (12) 

where a = —123.6 MeV, b = 70.4 MeV and cr = 2 corresponding to the stiff nuclear 
EOS of isospin-symmetric nuclear matter. The asymmetry pressure Pasy due to the isospin- 
dependent potential is 

Pasy = [SoiPo) ■ Po " 7^^^' " 8.5^^/=^] . (13) 

Shown in Fig. 4 is a comparison of the density derivatives of the three partial pressures. As 
an illustration we have selected the parameter set of 5 = 0.2, T = 5 MeV, 7=0.5 and 2. In 
the whole density range reachable in heavy-ion collisions at intermediate energies, the gradi- 
ent of the total pressure is overwhelmed by the contributions from the isospin-independent 
interaction and the kinetic pressure. The contribution from the isospin-dependent inter- 
action is relatively small. Nevertheless, it is still useful to note that both the symmetry 
pressure and its density gradient increase with the 7 parameter, especially at high densities. 
We thus expect to see a stronger flow signal with a larger 7 parameter, but the effect should 
be rather small. This expectation is verified by performing the standard analysis of the 
transverse collective flow ED 

< pJA > (y) = PiM. (14) 

where A{y) is the number of particles at rapidity y and is i^^ particle's transverse mo- 
mentum in the reaction plane. As shown in Fig. 5., a slightly stronger transverse flow signal 
is obtained with 7=2 than 7=0.5 as expected. Given the difficulties and uncertainties in 
determining the reaction plane, the observed effect of the isospin dependence of the nuclear 
EOS is too small to be practically very useful. To explore the isospin-dependent EOS, one 
thus has to look for more delicate observables, such as the neutron-proton differential flow 



as proposed in ref. [|13 



As a complementary approach which also has the advantage of not having to measure the 
collective flow for neutrons, we investigate the proton differential elliptic flow at midrapidity 
as a function of transverse momentum ||?1|J33 



1 N ^2 _ 2 

<-M>=j,i:-'^. (15) 

i=l Pit 

where is the total number of unbound protons in the rapidity range of —0.5 < 
{y/Uhearnjcms < 0.5. In this rapidity range the isospin fractionation is the strongest, i.e., 
the value of N depends most sensitively on the parameter 7, as shown in Fig. 3. The piy 
is the ith particle's transverse momentum perpendicular to the reaction plane. The value 
of V2{pt) is thus a differential measure of the second-order flow effect as a function of trans- 
verse momentum. The sign and magnitude of f 2 reflects the result of a competition between 
the early "squeeze-out" {piy dominates) perpendicular to the reaction plane and the later 
in-plane transverse flow of nucleons. The differential elliptic flow, especially at high pt, 
is expected to be more sensitive to the isospin dependence of the nuclear EOS than the 
flrst-order transverse flow. This is because all three partial pressures lead approximately to 
a similar difference 6ply = pi — p^ although their respective contributions to the value of 
< Px > OT < Py > is very different. Moreover, the early pressure created in the participant 
region is revealed more clearly by the value of V2{pt) at high transverse momenta. This is due 
to the fact that high pt particles can only be produced through the most violent collisions 
in the early stage of the reaction. This fact is almost universal in heavy-ion collisions at all 
energies, see, e.g., [^3|-|4^. These particles can only retain their high transverse momenta 



by escaping from the reaction zone along the direction perpendicular to the reaction plane 
without suffering much rescatterings. 

Shown in Fig. 6 are the predicted proton differential elliptic flow as a function of trans- 
verse momentum pt. The results presented here are obtained by using 20,000 Sn+Sn events 
in each case. Below the Fermi momentum of about 0.3 GeV/c the elliptic flow increases with 
Pt as predicted by the nuclear hydrodynamics P7| , ^ . This is because the low pt particles 



have undergone a sufficiently large number of rescatterings to reach the hydrodynamical 
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limit and to obtain an in- plane flow indicated by the positive value of V2{pt)- Particles with 
higher transverse momenta are mainly those emitted perpendicular to the reaction plane 
with negative values of V2. These particles carry more undisturbed information about the 
initial high density phase of the reaction. This effect becomes stronger gradually with the 
increasing transverse momentum. Therefore the v^ipt) value decreases with the increas- 
ing pt above the Fermi momentum. Indeed, it is seen that the elliptic flow at high pt is 
much more sensitive to the parameter 7 than the flrst-order transverse flow. Without the 
Coulomb repulsion the shadowing effect of the spectator nucleons lasts longer, it leads to 
a more stronger "squeez-out" and thus a more negative value of elliptic flow, especially at 
high Pt- As the pressure gradient increases by increasing the parameter 7, the in-plane flow 
becomes dominant over the "squeeze-out" in the whole range of transverse momentum. The 
value of V2{pt) thus increases with the parameter 7, particular in the high pt region. The 
sensitivity observed here is much stronger than that observed in analyzing the flrst-order 
transverse flow. Although very high pt particles are rather rare, it is practically relatively 
easy for nucleons to obtain a transverse momentum upto 2 times the Fermi momentum in 
heavy-ion collisions at intermediate energies. In the range of pt < 0.5 GeV/c, our results 
here indicate an unambiguously strong signal of the isospin dependence of the nuclear EOS 
far beyond the statistical error bars. 

In summary, within an isospin-dependent transport model for nuclear reactions involving 
neutron-rich nuclei, we analyzed the flrst-order direct flow and the second-order differen- 
tial eUiptic flow for free protons in the reaction of ^^^S'n +^^'^ Sn at a beam energy of 50 
MeV/nucleon and an impact parameter of 5 fm. The emphasis was placed on searching for 
observable signals of the isospin dependence of the nuclear EOS. We found that the differen- 
tial elliptic flow, especiaUy at high transverse of momenta, for midrapidity protons is much 
more sensitive to the isospin dependence of the nuclear EOS than the flrst-order direct flow. 
This flnding provides a useful guide for planning future experiments to study the isospin 
dependence of the nuclear EOS by using neutron-rich heavy-ion collisions at intermediate 
energies. 
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FIG. 1. Symmetry potential for neutrons (upper branches) and protons (lower branches) as a 
function of density for an isospin asymmetry of 5 = 0.2 and the 7 parameter of 0.5, 1.0 and 2.0, 
respectively. 
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FIG. 2. A scatter plot of the nucleon density of each cell versus its isospin asymmetry 6. 
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FIG. 4. Derivatives of the three partial pressures as a function of density with the isospin 
asymmetry 6 = 0.2 and the 7 parameter of 0.5 and 2. The sohd (dashed and long-dashed) line is 
for the isospin-independent (dependent) interaction pressure, while the dotted line is for the kinetic 
pressure at a temperature of 5 MeV. 



18 




FIG. 5. The average in-plane transverse momentum of free protons as a function of rapidity 
for the Sn+Sn reaction. 
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